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Q> ■ Abstract 



By exploiting the analyticity properties of the two-point current-current cor- 
relator we obtain numerical predictions for the e^e~ moments in terms of the r 
decay rate. We perform a partial resummation of the pertinent perturbative series 
, expansion by solving the renormalization group equation for Adler's function. Our 

predictions are renormalization scheme independent but depend on the order of the 
perturbative /3-function expansion. The analysis involves the unknown five-loop 
coefficient for which we give some new estimates. 



The interpretation of experimental strong interaction data in terms of QCD and the 
Standard Model requires higher and higher accuracy in the theoretical input There 
is a need to improve on the precision of the predictions of perturbative QCD by developing 
a framework that allows one to go beyond finite orders of perturbation theory. This is 
indispensable if one wants to make progress in precision fits to data of existing and 
especially future experiments. New high order results in fixed order perturbation theory 
are becoming more difficult to come by because of overwhelming technical difficulties. [] 
What is needed is to develop techniques that allow one to go beyond the existing finite 
order perturbative results. Attempts in this direction include finite order predictions based 
on the Fade approximation different optimizations of perturbation theory [§, and 
infinite resummation procedures based on particular properties of perturbation series like 
renormalon methods H. Another approach consists in the use of the renormalization 
scheme freedom to parameterize infrared contributions to physical observables in the 
renormalon approximation ||^, |10|, [ll|]. Latter approaches can serve as an alternative to 



the renormalon calculus. 

The analysis of the moments of the e~^e~ rate and, in conjunction with it, the analysis of 
the T decay rate has a long-standing history both in theory and experiment. The accuracy 
in the determination of the r decay rate and its decay characteristics has been continuously 
improving while there is now hope that there will be more precise data on the e~^e~ rate 



in the low energy domain in the near future [|T2|. Because of the availability of a large 
number of terms in the perturbative QCD expansion and the simplicity of the analyticity 
structure of the underlying Green's function, the analysis of the e^e~ annihilation process 
has advanced to a highly sophisticated stage. The analyticity structure is simple because 
the process is related to the two-point correlator of gauge invariant current correlators, 
the analyticity structure of which is determined by the Kallen-Lehmann representation. 

In a recent paper we advocated the idea to directly compare physical observables 
within fixed order perturbative QCD. Using this approach we determined moments of the 
e~^e~ annihilation rate in terms of the r decay rate |T^. This eliminates the problem 
of scheme dependence within finite order perturbative QCD (especially if the moments 
are taken at the same or comparable scale). In a second paper we exploited analyticity 
properties of the two-point current correlator to partially resum the perturbation series for 
a new analysis of 0] (see also In the present paper we combine the approaches 



of [|I^ and [|I^ and present the results of an analysis which expresses the moments of 
e~^e~ rate functions in terms of the r decay rate including resummation effects. We also 
discuss some general features of the solution of the renormalization group equation in the 
complex plane and speculate on estimates of higher order coefficients of the perturbation 
series in a renormalization group invariant manner. 

We closely follow the notation introduced in [|1^. Let us begin by defining moments 
of the e+e~ annihilation rate in terms of the spectral density R{s) according to 

/■So fl a / Q \^ 

R4so) = in + l) - - Ris). (1) 

Jo So VSq/ 

For convenience we have normalized the moments such that Rn{so) = 1 for R{s) = 1. We 
define reduced moment functions by factoring out moments of the partonic Born term 



^One of the recent achievements was the calculation of the four- loop /3-function and the four- loop 
anomalous mass dimension [|[. 



2 



contribution R^{so). The moments of the Born term contribution are given by the 
leading order in the strong interaction. Nonperturbative (power suppressed) corrections 
and corrections due to other interactions (i.e. electro-weak corrections in the case of the 
semileptonic r decay) may be absorbed in |]T6[. After factorization of the Born term 
moments the reduced moment functions r„(so) are given by 

i?„(so)=i?°(so)(l + ^r„(so)). (2) 

Our aim is to establish relations between different sets of observables (and not relations 
between observables and powers of the strong coupling constant). Further the relations 
should be independent of the choice of a particular renormalization scheme. It is therefore 
convenient to deal with scheme independent quantities from the outset. 

Similar to Eq. (0) we define the reduced Adler's function d^Q"^) 

D{Q') = A fl + id{Q')] (3) 
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which we take to be an effective coupling constant. In the Euclidean domain the reduced 
Adler's function d{Q^) has the expansion 

diQ^) = a(g2) + k,a\Q^) + ha^Q^) + ha\Q^) + ... (4) 

where a = 9as/4:7i. 

The running of the effective coupling constant d{Q'^) is determined by the renormal- 
ization group equation which we write as 

Q'^^(Q')=/?W)), (5) 

where the perturbative series expansion of the /3-function reads 

p{d) = -d'^il + pirf + p2d^ + p^d"" + ...). (6) 

Due to the fact that we have chosen the reduced Adler's function as our expansion param- 
eter, the coefficients are renormalization scheme independent quantities (see also [p!8|). 
Numerically they are given by [N^ = 3) 

64 

= _ ^ 0.790 
^ 81 

16531 728 , , , ,y 200 , , 

= ,37096148 4820288 12352 , _ 

59049 6561 ' 81 ^ ' 

-^C(5) + ^C(3)C(5) + 2h ^ 2h - 2.97953 (7) 
where we have taken the numerical values for the known coefficients in the MS scheme 



from [|I9[ together with the recently calculated four-loop /3-function coefficient given in P] . 
It should be kept in mind though that the final values of the coefficients are scheme 
independent even if they have been calculated in a specific scheme. The five-loop contri- 
bution ^3 in the coefficient ps is not yet known. Later on we shall present some estimates 
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for ks- We have set = 3 in the present apphcation. Here we are interested in moments 
of e~^e~ annihilation spectral density at the scale = m^. The extension of our approach 
to other scales is straightforward. 

The reduced Adler's function d{m'^) itself can be determined from experimental data 
in terms of the integral representation 



d{ml) = m. 



r{s)ds 

{s + miy 



where r{s) is the reduced spectral density. Note that after factorization of the Born term 
contribution, r(s) need not be positive definite anymore. One therefore may encounter 
dramatic cancellations in the evaluation of the integral in Eq. (^) at the cost of the 
precision with which d{rn^) can be determined. This is the reason why we do not use the 
integral respresentation Eq. (§) for our numerical estimates of d{m^). 

The standard technique of contour integration in the complex plane allows one to 
obtain closed formulas for the reduced moments r„ in terms of d(m^) = d-r and the 
coefficients pi appearing in the perturbative expansion of the /5-function in Eq. (P). The 
moments are given by 

rn{m1) = — f x^p{mlx)dx, (9) 



27rz J|x|=i 

where p{z) is the reduced vacuum polarization function. It is related to the reduced 
Adler's function by 

dm = -Q'^^fe^- (10) 



After integrating Eq. by parts using Eq. (|TUp, the moments can be represented by 

r„(m^) = rcirc(m^) + A„(m^). (11) 

The first term on the right hand side of Eq. ([TT|) represents the surface term contribution. 
It does not depend on n and is, in fact, the renormalization group improved spectral 
density (the discontinuity across the cut at s = ± iO) which can be calculated from 
d{Q'^) = d{mle'"^) using 

d{my^)d(j). (12) 



(ml] 
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27r 

The second term in the partial integration is given by 



Animl) = / e'("+i)<^rf(m2e'<^)# (13) 



27T 

and depends on the order n of the moment under consideration. Latter contributions 
tend to be numerically suppressed because of the oscillatory factor e*("+^)'^. Considered 
as functions of the variable n, the quantities A„ are oscillating functions. At discrete 
values for n they can give quite irregular contributions to the moments. Numerically 
these contributions are suppressed but are not negligible. 

We mention that the reduced moments cannot be approximated with any precision 
by an asymptotic expansion in dr- For small enough values of d^- the corresponding series 



converge WM. However, if we take the numerical values for dr derived below, they lie 



outside the circle of convergence and so the terms in Eq. ([TT|) cannot be represented as 
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series in dr- That is why we prefer to work directly with the integral representation given 
by Eqs. (0) and (pTSD and not with a series expansion. 

To leading order in the /3-function the reduced moments can be presented in explicit 
form. With f3{d) = —d"^ one can integrate the renormalization group equation in the 
complex Q^-plane using the starting value dr and obtains d{m'^e^'f') = dr{l + idT-(f))~^. The 
moments can then be represented as 



2, 1 , ^ , (-1)" r e*("+i)<^# 



r„(m^) = - arctan(7rrf^) + — / ——3—-. (14) 

TT ZTT J-TT 1 + ld-r(p 

If more terms are included in the /3-function expansion, we write the general result in 
symbolic form as 

r„(m^) = r„(rf^; p) = r„(rf^; pi, p2, Ps, • • •)• (15) 

When expanded for small dr and a given order of the /3-function expansion, these moments 
reproduce the known results of fixed order perturbation theory. However, here we consider 
a different approach which is not based on a series expansion. The functions r„(m^) can be 
computed from Eq. (^ to any given order of the /3-function expansion without ever having 
to invoke a series expansion. When one goes beyond the leading order approximation 
of the /3-function, the renormalization group equation (^ can still be solved, but the 
integrations in Eqs. (|T2|) and (|T3|) become unwieldy. In these cases we proceed in numerical 
fashion. 

In order to make contact with the semileptonic r decay rate we take a particular 
combination of moments, namely 

rr = 2ro{ml) - 2r2(m^) + r3(m^). (16) 

Using the partially integrated representation Eq. (pAj), this results in 

Triml) = r,,,,{ml) + 2Ao{ml) - 2A2{ml) + Asiml) = r,i,,{ml) + Ar{ml). (17) 

This formula holds for perturbative QCD in the case of massless quarks where axial- vector 
and vector contributions are equal. 

As a next step we invert Eq. (0) numerically to obtain numerical values for dr by 
solving the differential equation (HD with dr as starting value. This is done by a systematic 
trial and error procedure using as many iterations as necessary. The whole procedure is 
repeated for successive orders in the /3-function expansion. Using the experimental value 



0.487 ±0.011 g, we obtain 



4°^ = 0.3431 ± 0.0069, rf^ = 0.3466 ± 0.0081, = 0.3561 ± 0.0097 (18) 

for increasing orders in the /3-function expansion, where the superscript labels the order 
of the expansion. Once these values are known, we can predict the contributions of rdrc 
and A„ to the different reduced moments r„. Predictions for the different contributions 
using the leading, first and second order accuracy for the /3-function are given in the first 
three columns of Table |l|, whereas in Table ^ we present the results in terms of the full 
moments. 

For the orders n = 0, 1 and 2 of the moments the values in Table ^ agree with estimates 
given in |l^ on the basis of fixed order expansion in perturbation theory. We cannot give a 



prediction for the moment r_i/2 within the present approach because the weight function 
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3 




0.3431 


0.3466 


0.3561 


0.3667 


' CllC \ ' ' "7- J 


0.2619 


0.2396 


0.2352 


0.2354 


AoK) 


0.1255 


0.1369 


0.1411 


0.1448 


Ai(m2) 


0.0290 


0.0174 


0.0108 


0.0073 


A2(m2) 


0.0204 


0.0196 


0.0210 


0.0227 


AsK) 


0.0148 


0.0129 


0.0118 


0.0110 


A4(m2) 


0.0116 


0.0105 


0.0103 


0.0106 


A,(m2) 


0.2251 


0.2474 


0.2519 


0.2552 



Table 1: Moment contributions for increasing /^-function accuracy i 
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2.3444 


2.3347 


2.3344 


2.3380 
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2.2586 


2.2284 


2.2187 


2.2157 


R2 
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2.2509 


2.2304 


2.2277 


2.2295 


R3 


[ml) 


2.2459 


2.2245 


2.2195 


2.2190 


Ri 


[ml) 


2.2432 


2.2223 


2.2182 


2.2187 



Table 2: Full moments for increasing /3-function accuracy i 



used for constructing moments has to be analytically continued and must be single- valued 
in the vicinity of the cut along the positive semi-axis. This is the only singularity in the 
complex plane allowed for the vacuum polarization two-point function as given by the 
spectral representation. Though ^/z is an analytical function in the complex plane with 
a cut, its values on different sides of the cut have opposite signs. Because of this, one 
cannot extract the discontinuity of the vacuum polarization function across the cut. 

As a note we want to illustrate the advantage of the resummation technique in compar- 
ison with ordinary perturbation theory. This illustration is done by giving a very simple 
example concerning the really non-trivial question whether results can be obtained by 
directly keeping larger number of terms. Consider two observables given by perturbative 
series in the some given scheme, namely 

/(a) = a(l -a + a^ -...) = (19) 

1 + a 

and 

g{a) = a(l - 2a + 4a2 - . . . ) = (20) 
The functions /(a) and g{a) can be seen to be related by 

^(/) = ^ = /(l-/ + /'----)- (21) 

If we fit the right hand side of Eq. ([I9|) to an experimental value of about / = 0.6, we get 
a = 1.5. But for this value of the coupling, the series in Eq. ( |T^ ) diverges. So we cannot 
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get a from it without a proper resummation procedure that in this case is trivially given 
by the appended exact formula. Consequently we cannot get a prediction for g using the 
series in Eq. ( ^OD in terms of a. On the other hand, the direct relation in terms of the 
series in Eq. ( pT| ) converges perfectly and gives an unambiguous result for g in terms of 
measured /. Of course, for such an improvement to occur one has to analyze in detail 
the underlying theory and the origin of the series. Within our resummation procedure we 
are able to do so though we still could not answer whether one can reexpand moments 
through Vj. directly to get a convergent series. We could not prove the opposite either. 

Comparing Eq. ([ToD with Eq. (0) and Eq. (^OD with Eq. ([T5|), the example demon- 
strates that we can obtain a dependence between the reduced e'^e~ moments r„ and the 
r decay rate Vr without using an expansion in dr. This implicit dependence is shown in 
Fig. 1 for different moments. 

In contrast to our previous study [|13| the e~^e~ moments can be computed without hav- 
ing to perform a sophisticated analysis of divergent series and estimating the truncation 
errors. The only question that remains is the question of errors in the present approach. 
The statistical error resulting from the uncertainty in the experimental number for r^- can 
easily be taken into account, while for the perturbation series itself we suggest to take the 
difference between results for different /3-function accuracy as the resulting error. When 
analyzed along the contour in the complex Q^-plane, the worst pattern of convergence 
for the /3-function expansion is given in the vicinity of the Euclidean point = and 
reads 

/?(0.36) = -(0.36)^(1 + 0.79(0.36) + 1.035(0.36)2 + P3(0.36)^) 

= -(0.36)2(1 + 0.284 + 0.134 + 0.0467p3)- (22) 

To estimate the error of this expansion we have to estimate the value of ps. This is also 
necessary to obtain a feeling for the accuracy of the perturbative approximation for the 
/9-function in Eq. (^2|) . In the MS scheme, this coefficient is given by the coefficients of the 
/9-function up to the known coefficient P3 and the coefficients of the Adler's function 
up to the yet unknown five-loop coefficient k^. In the MS scheme there exist estimates 
for this quantity. They are essentially Fade estimates valid within this particular scheme 
and result in a value of = 2.17 or p3 = 1.36. Other estimates are based on various 
optimization procedures for the perturbation series. They give values close to the Fade 
estimate 0. 

Let us add to the above estimates of ps and present our own analysis. The first estimate 
is a Fade approximation for the /5-function itself, which gives pa = p|/pi = 1.3548. We 
obtain a second estimate by considering a one parameter subgroup of the renormalization 
group which leaves the /3-function invariant. It is given by 

d'iml) = d{e^ml) = d^ - id^ + (7^ - P17) - (7' - \pil' + P27X + • • • (23) 

and thus expresses ps as a function of 7. Although the overall value of pa is scheme inde- 
pendent by definition, one introduces a scheme dependence for pa through the estimation 
procedure because k^ is scheme dependent. The dependence of ps in terms 7 is shown in 
Fig. 2. We see that the value of ps = 1.36 in the MS scheme (i.e. for 7 = 0) is not stable 
against small variations of 7. There is, however, a region where ps is almost independent 
of 7 and yet close to the value given in the MS scheme. Choosing a scheme in the stability 
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region is known as the principle of minimal sensitivity (PMS) |2T[ and works well in a 



number of applications (see for instance p2|). Using this principle, we obtain the value 



P3 = 2.4530 in a wide range around the value 7 = —1.3288. It happens that this choice 
is close to the so-called G-scheme with 7 = —2, where p3 = 2.0518. Although the 
above dependence is not the most general variation within the renormalization group, we 
believe that it gives an additional support for the obtained value of p^. 

The last estimate may look a bit extravagant, nevertheless the result is consistent with 
that of the previous ones. We fix a scheme such that the first few terms of the reduced 
Adler's function d^Q"^) are represented by a pure geometric series 

diQ') = aGs(Q')(l + kacsiQ') + k^^lsiQ^) + fc'«Gs(Q') • • • ) (24) 

with coefficients of the /5-function given in the MS scheme. This is always possible up 
to order k"^ because of the freedom of the one-dimensional reparametrization invariance 
given by Eq. (^). Numerically one finds 

c?(Q') = aGs(Q')(l-0.1917aGs(Q')+0.0367a2.g(Q2) + (A;3_2.602)a^s(<9') + --- ) (25) 

where fcs is the unknown five-loop coefficient given in the MS scheme. By demanding 
a geometric series behaviour for Eq. ( p5|) one calculates = 2.595. As an average of 
our three estimates we quote ps = 2.0 it 0.5. With this final estimate for we obtain 
dr = 0.3667 ± 0.0120 and the results given in the last columns of Table and Table 

The statistical error due to the input uncertainties in is about 3% for all moments. 
This error is larger than the change resulting from adding the next (estimated) term of 
the /3-function expansion. We therefore conclude that up to this order in the /3-function 
the perturbative expansion for the /3-function does not seem to limit the accuracy of 
the resummed predictions for the e~^e~ moments. The main uncertainty comes from the 
experimental error in the semileptonic r decay rate. 

To conclude, we have obtained numerical predictions for the moments of the e~^e~ 
annihilation rate in terms of the known r decay rate. We exploited the analyticity prop- 
erties of the two-point current-current correlator to perform a partial resummation of 
the perturbative series relating the two sets of observables. Our predictions for the e~^e~ 
moments are renormalization scheme independent but depend on the order of the pertur- 
bative /3-function expansion. We have attempted to estimate the error in our prediction 
for the e~^e~ moments by estimating an unknown five-loop piece in the /5-function expan- 
sion. Using this estimate, we found that the error in our predictions is dominated by the 
experimental error of the r decay rate. 
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Figure Captions 

Fig. 1: Implicit dependence of the reduced e"'"e~ moments ro, ri, r2 and on the 
semileptonic r decay rate 

Fig. 2: Dependence of the /3-f unction coefficient ps on the subgroup parameter 7 
which specifies the choice of the renormalization scheme 
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